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The variety of isotopes in cosmic rays allows us to study different aspects of the processes that cosmic rays 
undergo between the time they are produced and the time of their arrival in the heliosphere. In this paper we 
present measurements of the isotopic ratios ^H/*He, ^He/'^He, ^Li/^Li, ^Be/(^Be+^"Be) and ^"B/^^B in the 
range 0.2 — 1.4 GeV of kinetic energy per nucleon. The measurements are based on the data collected by the 
Alpha Magnetic Spectrometer, AMS-01, during the STS-91 flight in 1998 June. 

Subject headings: acceleration of particles — cosmic rays — nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

Cosmic rays (CRs) detected with kinetic energies in the 
range from MeV to TeV per nucleon are believed to be pro- 
duced by galactic sources. Observations of X-ray and 7-ray 
emission from galactic sites such as supernova remnants, pul- 
sars or stellar winds reveal the presence of energetic particle 
acceleration mechanisms occurring in such objects. The sub- 
sequent destruction of these accelerated nuclei {e.g., p. He, 
C, N, O, Fe) in the interstellar medium gives rise to sec- 
ondary species that are rare in the cosmic ray sources, such as 
Li, Be, B, sub-Fe elements, deuterons, antiprotons, positrons 
and high energy photons. The relation between secondary 
CRs and their primary progenitors allows the determination 
of propagation parameters, such as the diffusion coefficient 
and the size of the diffusion region. For a recent review, see 
Strong et al. (2007). 

Along with the ratios B/C and sub-Fe/Fe, it is of great im- 
portance to determine the propagation history of the lighter 
H, He, Li and Be isotopes. Since ■^H and ^^He CRs are mainly 
produced by the breakup of the primary *He in the galaxy, 
the ratios ^H/*He and "^He/^He probe the propagation history 
of helium (Webber 1997). The isotopes of Li, Be and B, all 
of secondary origin, are also useful for a quantitative under- 
standing of CR propagation. The relative abundances and iso- 
topic composition of H, He, Li, Be and B, therefore, might 
help to distinguish between the propagation models and give 
constraints to their parameters (Moskalenko et al. 2003). 

Low energy data (<200 MeV nucleon"^) on CR isotopic 
composition come mainly from space experiments such as the 
HET telescopes on VOYAGER 1 and 2 (Webber et al. 2002), 
the Cosmic Ray Isotope Spectrometer (CRIS) on the Ad- 
vanced Composition Explorer (ACE) satellite (de Nolfo et al. 
2001), the ULYSSES high energy telescope (ConneU 
1998) and the HKH experiment on the International Sun- 
Earth Explorer (ISEE) spacecraft (Wiedenbeck & Greiner 
1980). Light nuclei data at higher energies (up to few 
GeV nucleon^ ^) have been measured by balloon borne mag- 
netic spectrometers including IMAX (Reimeretal. 1998), 
ISOMAX (Hamsetal. 2004), SMILI (Ahlenetal. 2000), 
BESS (Wang et al. 2002), Inteplanetary Monitoring Platform 
(IMP) experiment (Garcia-Munoz 1977) and the Goddard 
Space FHght Center (GSFC) bafloon (Hagen et al. 1977). 

AMS-01 observed CRs at an altitude of ^^380 km during 
a period, 1998 June, of relatively quiet solar activity. It col- 
lected data free from atmospheric induced background. In this 
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paper we present measurements of the ^He/'*He ratio over the 
kinetic energy range 0.2 — 1.4 GeV per nucleon, and the aver- 
age values of the ratios ^LifU, ■^Be/(^BeH-i°Be) and i"B/"B 
over the same energy range. The ratio ^H/^He is also pre- 
sented. 

2. THE ALPHA MAGNETIC SPECTROMETER 

The Alpha Magnetic Spectrometer (AMS) is a particle 
physics instrument designed for the high precision and long 
duration measurement of CRs in space. The AMS-01 precur- 
sor experiment operated successfully during a 10 day flight on 
the space shuttle Discovery (STS-91). 

The spectrometer was composed of a cylindrical permanent 
magnet, a silicon micro-strip tracker, time-of-flight (TOF) 
scintillator planes, an aerogel Cerenkov counter and anti- 
coincidence counters. The performance of AMS-01 is de- 
scribed elsewhere (Aguilar et al. 2002). 

Data collection started on 1998 June 3. The orbital incli- 
nation was 51°. 7 and the geodetic altitude ranged from 320 
to 390 km. The data were collected in four phases: (a) 1 day 
of check out before docking with the MIR space station, (b) 
4 days while docked to MIR, (c) 3.5 days with AMS point- 
ing directions within 0°, 20° and 45° of the zenith and (d) 0.5 
days before descending, pointing toward the nadir. 

The acceptance criterion of the trigger logic in the AMS- 
01 instrument was a four-fold coincidence between the sig- 
nals from the four TOF planes. Only particles traversing 
the silicon tracker were accepted. Events crossing the anti- 
coincidence counters or producing multiple hits in the TOF 
layers were rejected. A prescaled subsample of 1 out of 1000 
events was recorded with a dedicated minimum-bias configu- 
ration. This "unbiased trigger" required only the TOF coinci- 
dence. 

The AMS-01 mission provided results on cosmic 
ray protons, helium, electrons, positrons and light nu- 
clei (Aguilar et al. 2002, 2007, 2010). During the flight, a 
total of 99 million triggers were recorded by the spectrometer, 
with 2.85M helium nuclei and nearly 200,000 nuclei with 
charge Z > 2. 

3. DATA ANALYSIS 

The identification of cosmic ray nuclei with AMS-01 was 
performed through the combination of independent measure- 
ments provided by the various detectors. The particle rigid- 
ity, R, (momentum per unit charge, pc/Ze) was provided by 
the deflection of the reconstructed particle trajectory in the 
magnetic field. The velocity, (3 = v/c, was measured from 
the particle transit time between the four TOF planes along 
the track length. The reconstruction algorithm provided, to- 
gether with the measured quantities R and /3, an estimation 
of their uncertainties 6R (from tracking) and 6/3 (from tim- 
ing), that reflected the quality of the spectrometer in per- 
forming such measurements. The particle charge magnitude 
\Z\ was obtained by the analysis of the multiple measure- 
ments of energy deposition in the four TOF scintillators up 
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io Z = 2 (Aguilar et al. 2002) and the six siHcon layers up 
to Z = 8 (Aguilar et al. 2010). The particle mass number. A, 
was therefore determined from the resulting charge, velocity 
and rigidity: 

A^^VT^ (1) 

where m„ is the nucleon mass. 

The response of the detector was simulated using the AMS 
simulation program, based on GE ANT -3. 21 (Brun et al. 
1987) and interfaced with the hadronic package RQMD (rel- 
ativistic quantum molecular dynamics (Sorge 1995)). The 
effects of energy loss, multiple scattering, nuclear interactions 
and decays were included, as well as detector efficiency and 
resolution. After the flight, the detector was extensively cali- 
brated at GSI, Darmstadt, with ion beams (He, C) and at the 
CERN-PS, Geneva, with proton beams. This ensured that the 
performance of the detector and the analysis procedure were 
thoroughly understood. 

Further details are found in Aguilar et al. (2002) and refer- 
ences therein. 

3.1. Helium Isotopes 

Given the large amount of statistics available for Z = 2 
data, we considered only the highest quality data collected 
during the post-docking phase (c) and only while pointing to- 
ward the zenith. Data taken while passing near the South At- 
lantic Anomaly (latitude: 5° - 45° S, longitude: 5° - 85° W) 
were excluded. Only events taken when the energy interval 
0.2 — 1.4 GeV nucleon^^ was above the geomagnetic cutoff 
for both the isotopes ■^He and "'He were kept; this corresponds 
in selecting the orbital regions with the highest geomagnetic 
latitudes, 0m, roughly |0m| ^ 0.9. 

Furthermore, the acceptance was restricted to particles 
traversing the detector top-down within 30° of the positive 
z-axis. Events with poorly reconstructed trajectories were 
rejected through quality cuts on the associated or con- 
sistency requirements between the two reconstructed half 
tracks (Aguilar et al. 2010). To avoid biasing the recon- 
structed mass distributions, no cuts on the consistency of the 
TOF velocity versus tracker rigidity measurements were ap- 
plied. We required that the velocity was measured with hits 
from at least three out of four TOF planes and that the rigidity 
was reconstructed with at least five out of six tracker layers. 
Approximately 18,000 nuclei with charge Z ~ 2 were se- 
lected in the energy range 0.2 — 1.4 GeV nucleon"^. The 
charge was determined from the energy depositions in both 
the TOF and tracker layers. The kinetic energy per nucleon 
was measured with the TOF system, i.e., through the velocity 
p. In the considered energy range, the TOF energy resolution 
is comparable to that of the tracker 

The selected data are shown in Fig. 1 distributed in the 
R) plane. The two curves represent the exact relation be- 
tween velocity /3 and rigidity R for a Z = 2 nucleus of mass 
number A = 3 (dashed line) and A = 4 (solid Une), which is: 



/3 



mn(_ 
ZeR 



-1/2 



(2) 



The large dispersion of the measured data, apparent from 
Fig. 1, indicates a relatively poor mass resolution in the sepa- 
ration of the two mass numbers. Under these conditions, any 
event-to-event separation {e.g., through a mass cut) is clearly 
inapplicable. In addition, the distribution of the reconstructed 
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Fig. 1. — Distribution of the measured velocity, /3, as a function of tlie 
reconstructed rigidity, R, for Z = 2 nuclei. The two lines represent the exact 
relationship of /3 to i? for the two isotopes ^He (dashed line) and *He (solid 
line). 



mass numbers (Eq. 1) exhibits asymmetric tails; so the stan- 
dard Gaussian fit method (Seo et al. 1997) is not appropriate 
for describing the observed mass response of the instrument. 

In order to determine the isotopic ratios, it was therefore 
necessary to develop a comprehensive model for the complete 
response of the instrument to different masses. The mass res- 
olution is influenced by the intrinsic time resolution of the 
TOF system and by the bending power of the magnet coupled 
with the intrinsic spatial resolution of the tracker. Physical 
processes such as multiple scattering, energy losses and in- 
teractions along the particle path also contribute in shaping 
the reconstructed mass distributions. Thus, we modeled the 
AMS-01 mass response by means of our Monte Carlo (MC) 
simulation program, which includes all the aforementioned 
physical effects as well as the instrumental readout, providing 
a realistic description of particle tracking and timing on an 
event-by-event basis. The program was also tuned with data 
collected during the test beams. The resulting rigidity resolu- 
tion 5R/ R and velocity resolution SP/ [3 are shown in Fig. 2 
for test beam data (filled circles) and MC events (histograms). 
It can be seen that the mass resolution, approximately given 
by: 

-a) = -j + UJ ' 

was correctly simulated as the MC agrees with the data within 
-2%. 

Using a sufficiently large number of simulated events of 
■^He and ''He and with the ratio ■^He/*He of the detected events 
as a free parameter, we determined the best composition fit be- 
tween the simulated mass distributions and the measured one. 
In these fits the overall normalization, M, was also a free pa- 
rameter In principle, M should be fixed by the data, namely 
by the number of entries, Ae, of each mass histogram. De- 
viations of TV from its expected value A/e may indicate the 
presence of an unaccounted background, e.g., from charge 
misidentification. 

The results of this procedure are shown in Fig. 3, where 
the agreement between the measured mass histograms (filled 
circles) and the simulated histograms (fines) turned out to be 
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Fig. 2. — a) Rigidity and b) velocity resolutions of ttie AMS-01 tracker 
and TOF estimated with measured data from the test beam with E = 2 
GeV nucleon^^ helium nuclei. Data (circles) are compared with the MC 
simulation (histograms). The MC entries are normaUzed to the data. 



TABLE 1 

Fit results and correction factors for the ratio 

^HE/'*HE between 0.2 AND 1 .4 GEV OF KINETIC ENERGY PER 
NUCLEON. 



Energy 


Events 


Fit Results 


X^/df 


<5A/A 


ACorr 


FCorr 


0.20-0.30 


2,660 


0.125 ±0.011 


32.3/29 


13.1% 


1.12 


0.97 


0.30-0.44 


3,553 


0.158 ±0.096 


51.1/29 


12.2% 


1.05 


0.98 


0.44-0.64 


3,867 


0.182 ± 0.094 


65.8/34 


11.8% 


1.00 


0.98 


0.64-0.95 


4,142 


0.211 ± 0.098 


62.7/30 


12.2% 


0.99 


0.98 


0.95-1.40 


3,813 


0.223 ± 0.012 


55.0/33 


13.9% 


0.99 


0.98 



very satisfactory. The fits of the '^He/^He mass composition 
ratios gave unique minima in all the considered energy bins. 
The uncertainties associated to these ratios were directly de- 
termined from the l-tr uncertainties in the statistics of the 
fitting procedure. The fitting method was cross checked 
with the Maximum Likelihood method. The two methods 
gave the same results and very similar uncertainties. Dou- 
ble Gaussian fits were also performed in order to provide the 
corresponding mass resolution, (5 A/A, for each energy bin, de- 
fined as the ratio between the width and the mean'. The fitted 
'^He/^'He ratios for all the considered energy bins from 200 
MeV nucleon^'^ to 1.4 GeV nucleon^'^ are listed in Table 1, 
together with the x^/df values, the number of events, the mass 
resolution and correction factors discussed below. 

3.2. Top-Of-Instrument Corrections 

The measured mass distribution of Fig. 3 was fitted with 
an MC sample of mixed ■^He and ''^He that were sent through 
the same analysis chain (trigger, reconstruction and data se- 
lection) as the data. The free parameter is the ratio ■^He/^He 
of the two mass distributions corresponding to the recorded 
events. Hence, small corrections have to be performed in or- 
der to extract the ratio of interest, namely the ratio entering 
the instrument. These are referred to as Top-Of-Instrument 
(TOI) corrections. 

' Equal mass resolutions were obtained for "^He and *He within 0.1% at 
all energies. Table 1 provides the mean values. 



3.2.1. Acceptance Corrections 

The detector acceptance. A, includes trigger efficiency, re- 
construction efficiency and selection efficiency. The accep- 
tance was calculated using our Monte Carlo simulation pro- 
gram. Nucleus trajectories were simulated in the energy range 
^ 0.05 — 40 GeV nucleon^^. They were emitted downward 
from a square of length 3.9 m placed above the detector In 
total, 110 million ■^He and 550 million *He nuclei were sim- 
ulated. The physical processes involved and the detector re- 
sponse are very similar for the two isotopes, i.e., the resulting 
acceptances are quite similar in magnitude. The contributions 
from the detector acceptance mostly cancel in the ratio. The 
associated systematic errors also cancel. The only important 
factor in determining the isotopic ratios is the knowledge of 
any isotopic dependent effects in the detector response. Mass 
dependent features are expected from the following effects: 

• Rigidity threshold. The instrument acceptance is rigid- 
ity dependent, because the tracks of slower particles 
are more curved, and it is less likely that they pass 
through both the upper and lower TOF counters and 
the tracking volume. Since, at the same kinetic en- 
ergy per nucleon, the lighter isotope ■^He have lower 
rigidity than the heavier "^He, the resulting acceptance, 
particularly at lower energies, is lower for ■^He. Above 
0.2 GeV nucleon^'^, the rigidity threshold affects the 
ratio by less than ~ 1%. 

• Multiple scattering. Coulomb scattering is slightly 
more pronounced for lighter particles. Since multiple 
scattering affects the event reconstruction and selec- 
tion efficiency, the acceptance for the lighter isotope is 
smaller in the lowest energy region. The multiple scat- 
tering effect amounts to ^ 10% at E~0.2 GeV/n and 
decreases with energy, down to ^ 1 % at 1 GeV/n. 

• Nuclear interactions. The attenuation of cosmic rays 
after traversing the TOI material is isotope dependent 
and closely related to the inelastic cross sections, a^^^, 
for the interactions in the various layers of the detec- 
tor material. For -^He and ^He, the attenuation due to 
interactions differs by ^ 2% or less. 

Fig. 4 shows the ratio of the two acceptances for MC events 
after the successive application of the trigger, reconstruction 
and then selection cuts. Deviations are appreciable below 0.4 
GeV nucleon^ and mainly due to the event selection (filled 
squares). This indicates the dominance of the "multiple scat- 
tering effect" mentioned above, because the selection cuts 
acted against events with large scattering angles. 

These mass dependent features, due to the particle dynam- 
ics in the detector, did not appreciably influence the AMS-01 
trigger system. In Fig. 5 we report the reconstructed mass dis- 
tribution using data collected with the unbiased trigger (§2). 
The comparison of such unbiased data (circles) with a\\ Z ~2 
data collected from the flight (solid Une, normalized to the 
unbiased data entries) shows no significant difference in the 
mass distribution. 

3.2.2. Nuclear Interactions 

Table 2 lists the material between the top of the payload 
and the tracker: a Multi-Layer Insulation (MLI) blanket, a 
Low Energy Particle (LEP) shield and two TOF layers of 
plastic scintillators supported by a honeycomb structure. In 
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Fig. 3. — Reconstructed mass distributions of Zrec = 2 events from flight data (solid circles) and from MC generated isotopes of '^He (long-dashed lines), 
*He (short-dashed lines) and their sum (solid lines). Distributions are shown in five energy intervals from a) to e) and over the entire range /) between 0.2 and 
1.4GeVnucleon~^. 
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total there were ^ 5 g/cm^ of material above the tracking vol- 
ume. The RQMD interface used in the AMS sim- 
ulation program provided a simulation of all the high energy 



TABLE 2 

Material above the tracker. The 
column density is averaged over the 
angle of incidence. 



Detector Element 


Composition 


Amount 


MLI theiTnal blanket 


Cs H4 O2 


0.7 g/cm? 


LEP shield 


C 


1.3 gjcm? 


TOF scintillator 


C/H= 1 


2.1 gjcvp? 


TOF support structure 


Al 


1.0 gjcrn} 



hadronic collisions involving deuterons, ■^He, "'He and heav- 
ier ions. These effects give an appreciable contribution to the 
total acceptance of §3.2.1. The survival probability of •^He 
(*He) after traversing the TOl material of Table 2 varies be- 
tween -90% (-88%) at -0.2 GeV nucleon^i and -86% 
(- 86%) at - 1.4 GeV nucleon^i. 

More critical is the fragmentation of ^He into '^He, that re- 
quired a dedicated correction. In this process, if only a neu- 
tron is "stripped" above the tracker, the event is recorded as 
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Fig. 5. — a) Reconstructed mass distribution of Z = 2 events collected 
with the unbiased trigger (filled circles) in comparison with the Z = 2 mass 
histogram obtained with all the data (histogram), b) Ratio of the two his- 
tograms (squares); the horizontal dashed line is shown for reference. 



a clean "^He event. The measured "^He/^He ratio is then dis- 
torted by incoming "'He that spill over into the ■^He mass dis- 
tribution. Note that the simulated mass distributions of Fig. 3 
are referred to the particle identities within the tracking vol- 
ume, i.e., the ■^He mass histograms of the figures (long-dashed 
lines) also contain the "extra" "^He nuclei generated as frag- 
mentation products of ^He. Assuming that the kinetic energy 
per nucleon is maintained in the mass changing process, the 
ratio has been corrected for this effect. For each considered 
energy interval, the ratio jy between the "extra" "^He and the to- 
tal number of detected ^He was estimated. The isotopic ratio 
■^He/^He resulting from the composition fit Ai is then related 
to the TOI ratio TZ through the relation: 



M 



A3 4)3 + 04 

Ai (j>i — rjAi 04 



1 



1-77 



(4) 



where 03 and 04 are the incident (TOI) intensities of the two 
isotopes {TZ = 03/04) and A3 and Ai the corresponding ac- 
ceptances. Inverting Eq. 4 we obtain the TOI ratio: 



n 



Ai 
A3 



1 



77 



M 



M 



(5) 



We then define the TOI correction factors as 
ACorr = ^4/^3 for the acceptance, and FCorr = ^—V^m 
for the fragmentation. Note that ACorr is the quantity shown 
in Fig. 4 (filled squares). These values, to be applied as 
multiplicative factors to the fitted ratio, are also listed in 
Table 1. Such corrections are affected by <3% uncertainties 
in total, associated with the various physical and instrumental 
effects discussed here and in §3.2.1. All these errors and their 
role in the '^He/^He ratio are reviewed in §3.3. 

3.2.3. 5-Ray Emission 

The effect of (5-rays was included in our MC simulation. 
In our previous work (Aguilar et al. 2010), it was noted that 



energetic knock-on electrons affect the total acceptance at 
high energies. The production of (5-rays is proportional to 
the square of the primary particle charge, and the maximum 
energy of the produced (5-rays, E^^^, is proportional to the 
primary particle energy; for a nucleus of momentum Mj^c, 
approximately: 

Ei^^ = 2m,^p^^^ (6) 

For high energy nuclei {E > GeV nucleon^ '^), the emitted 
electrons can reach the anti-coincidence counters and veto the 
event, leading to an energy and charge dependent trigger effi- 
ciency. At lower energy, the (5-rays curl up inside the tracking 
volume, affecting the reconstruction efficiency. The influence 
of 5-rays below 1.4 GeV nucleon^ is negligible in this anal- 
ysis and their effect has no significant difference between iso- 
topes at the same energy, as M ^ m,,. 

3.2.4. Background 

The Z = 2 charge separation from Z < 2 and Z > 2 
samples was studied with MC simulations and inflight data 
of e~, p. He, and heavier ions. Proton (ion) beam data at 
CERN-PS (GSI) provided additional vaHdation at 0.75, 2.0, 
3.6 and 8 GeV nucleon^'^ of kinetic energy (Alcaraz et al. 
1999; Aguilar et al. 2002, 2010). The main potential source of 
background to the helium sample was protons and deuterons 
wrongly reconstructed as Z = 2 particles. Using the single 
TOE system or the single silicon tracker only, it can be seen, 
using flight data, that the probability of a Z = 1 particle to 
be reconstructed as Z = 2 is below 10^"^ thus affecting the 
helium sample of Using the combined measurements 

obtained from both the detectors, the probability of the wrong 
charge magnitude was estimated to be ^ 10^^ over all ener- 
gies (Alcaraz et al. 1999). Background from the less abundant 
Z > 2 particles is completely negligible compared to the sta- 
tistical uncertainties in the He sample (the He/Li ratio is ~ 200 
at the considered energies). 

3.2.5. Energy Losses and Resolution 

Charged cosmic rays that traverse the detector lose energy 
in the material above the tracking volume. The total energy 
loss has an appreciable effect on the lowest energy bins. The 
energy losses by Z = 2 nuclei were estimated and parameter- 
ized with the MC simulation program. We made an event-by- 
event correction to our data according to the average losses. 

Once the above corrections are performed, the relation be- 
tween the reconstructed energy of detected particles, Erec, and 
their true energy, Emi, is still affected by the finite resolution 
of the measurement. The probability of a bin-to-bin migra- 
tion P(i?REc|i?Toi) for He was estimated to effect only adjacent 
energy bins, to be symmetric and barely isotope dependent. 
Through these matrix elements, we estimated that the mea- 
sured ^He/^He ratio has an uncertainty of 1-3% due to the 
resolution. 

3.3. Uncertainty Estimate 

In Fig. 6, we summarize the various sources of uncertainty 
in the measurement of the '^He/"'He. Errors are organized in 
four categories: 

1. Mass fit. The dominant source of uncertainty (-^5- 
—9%) is that associated to fits on the mass distributions 
(§3.1). The errors were directly determined from the 
1-cr uncertainties in the statistics of the fitting pro- 
cedure. These errors are due to statistical fluctuations of 
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TABLE 3 

Uncertainty summary for the measured isotopic ratios. The various contributions are described in S3. 3. 



Error 



merHe vs Energy (GeV/n) 



Ratios in 0.2-1.4 GeV/n 



Type (Effect) 


0.2-0.3 


0.4-0.44 


0.44-0.64 


0.64-0.95 


0.95-1.4 


^He/^'He 


6Li/7Li 


7Be/9+"'Be 


it'B/iiB 


Mass Fit ((5 A/A & Statistics) 


8.9 % 


6.1 % 


5.2 % 


4.7 % 


5.3 % 


4.5 % 


9.0% 


15.9 % 


12.2 % 


Normalization {N-Ne) 


1.0% 


0.9 % 


1.3 % 


1.1 % 


1.9% 


1.5 % 


2.3 % 


3.0 % 


2.9%) 


Normalization (Acceptance) 


1.3% 


1.2 % 


1.1 % 


1.1 % 


1.0% 


0.5 % 


0.5 % 


0.6% 


0.5 %) 


Interactions (Inelastic) 


1.6% 


1.6 % 


1.8 % 


1.9% 


2.0% 


2.0% 


2.1 % 


2.7 % 


2.4%. 


Interactions (Fragmentation) 


2.8% 


2.2 % 


1.8 % 


2.5 % 


2.3 % 


2.7 % 


2.4% 


2.6% 


2.4 %£. 


Resolution (5/3//3) 


2.6% 


1.7 % 


1.4% 


0.8 % 


1.1 % 


0.6% 


0.6% 


0.6% 


0.6 %£. 


Total Uncertainty 


10.0 % 


7.0% 


6.2 % 


5.9% 


6.6% 


5.8 % 


9.8% 


16.7 % 


13.0 % 
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Fig. 6. — Relative eiTors on the '*He/''He ratio measurement as a function 
of the kinetic energy per nucleon. The total error (filled squares) is obtained 
by the sum in quadrature of all the other contributions. The lines are to guide 
the eye. 



the measured data and the inability of the spectrometer 
to separate the different masses within the mass resolu- 
tion SAJA. 

2. Normalization. As discussed in §3.1, two parameters 
A4 (ratio) and J\f (normalization) were fitted. An ideal 
fit should lead to M equal to the number of measured 
events, TVe- We took the relative difference M-Ne as 
a source of systematic error (^ 1%). Another contribu- 
tion (^ 1%) is due to the acceptance correction factors 
estimated with our MC simulation program (§3.2.1). 

3. Interactions. Our results rely partially on hadronic in- 
teraction models, as discussed in §3.2.2. Similarly to 
Wang et al. (2002), we assumed an uncertainty of 10% 
in the inelastic cross sections, which corresponds to 
~2% of systematic uncertainty. For the fragmenta- 
tion channel ^He^'^He, we assumed an uncertainty 
in the associated cross section equal to that cross sec- 
tion, obtaining an uncertainty of 2— 3% of the measured 
'^He/'*He ratio. Uncertainties in the material thickness 
were found to be negligible. 

4. Resolution. As discussed in §3.2.5, our measurement is 
affected by the finite energy resolution of the TOF sys- 
tem. A systematic uncertainty of 1—3% was estimated 
to account for this effect. 



The overall error (filled squares in Fig. 6) is taken to be the 
sum in quadrature of the different contributions. All these 
uncertainties are also reported in Table 3. 

3.4. Lithium, Beryllium and Boron 

In our previous work (Aguilar et al. 2010), the lithium iso- 
topic composition was determined between 2.5 and 6.3 GV of 
magnetic rigidity. Here we present a unified analysis of the 
lithium, beryllium and boron isotopes between 0.2 and 1.4 
GeV of kinetic energy per nucleon. In this measurement, we 
followed the same procedure as for the helium analysis. All 
the steps described in §3.1 were repeated for the study of the 
ratios ^Li/'^Li, ^Be/(9BeH-i°Be) and i°B/"B. In this section 
we outline the essential parts of the Z > 2 analysis. 

The capability of the spectrometer to separate isotopes 
close in mass was more critical for Z > 2 and the charge iden- 
tification capabilities were limited by the use of tracker infor- 
mation only. However, the most limiting factor for the Li-Be- 
B study was the statistics. Hence, we performed the measure- 
ment with just one energy bin from 0.2 to 1.4 GeV nucleon^ ^, 
and also included data from the MIR-docking (b) and post- 
docking, non-nadir pointing (c) phases (§2). As in our 
previous work, for data collected during phase (b), a geo- 
metric cut on the MIR shadow was applied to the accep- 
tance (Aguilar et al. 2010). The geomagnetic regions consid- 
ered and the event selection criteria were the same as for the 
helium analysis. Only four hits were required in the tracker, 
compared to five in the helium analysis. The particle charge 
was assigned using the identification algorithm described in 
Aguilar et al. (2010) and Tomassetti (2009), that was specifi- 
cally optimized for the Z > 2 species. 

Results from the mass composition fit are shown in Fig. 7. 
For comparison, the measurement was also performed on the 
average ratio of ■^He/^He over this energy range. The large 
statistical fluctuations of the Z > 2 data are apparent from the 
figure, in particular for the less abundant beryllium isotopes 
(400 events in total). However, the agreement between the 
measured masses (filled circles) and the simulated histograms 
(solid lines) was satisfactory. The TOI corrections to the mea- 
sured composition followed the procedure described in §3.2. 
In contrast to He, the Z > 2 acceptances were found to be 
a bit smaller for the heavier isotopes (^Li, '^'^"Be and ^^B) 
than the lighter ones (^Li, ^Be and ^"B), indicating the dom- 
inance of nuclear interactions over other effects (see §3.2.1). 
The Z > 2 mass resolutions, (5A/A, were found to be ^2,% 
larger than for the Z = 2 case, reflecting the slight charge de- 
pendence of the spectrometer performance in particle track- 
ing and timing. Corrections for fragmentation were also per- 
formed, considering the channels ^Li— >^Li, ^" '^Be— >^Be and 
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Fig. 7. — Reconstructed mass distributions of events reconstructed as a) lielium, b) lithium, c) beryllium and d) boron from data (filled circles) and MC samples 
of ^Li— '^Li, ''Be— '■'Be and ^"B— ^^B (lighter isotope: long-dashed lines; heavier isotope: short-dashed lines, sum: solid line). In c) the dotted line indicates a 
contribution of ^"Be, see the text. 
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Fig. 8. — Contour levels of the fitting paraineter ^Be/('-'Be-H^Be) and 
^"Be/^Be corresponding to la, 2a, and 3a bounds of uncertainty in the 
statistics. 

^^B— >^°B. A summary of fit results, TOI corrections and 
resolutions is given in Table 4. 

The errors were estimated as discussed in §3.3. The domi- 
nant error is that arising from the fit procedure: the small 
Z > 2 statistics and the broader mass distributions led to 
less constrained composition fits. Background from helium 
was estimated not to affect the lithium measurement, as for 
He— Li charge separation the TOF information was still us- 
able (Aguilar et al. 2010). More critical was the contamina- 



TABLE4 

Fit results and correction factors for the ratios 
^HE/*m, ^Li/^Li, ^Be/(Sbe-h10be) and iob/"B in the range 
0.2 — 1.4 GeV of kinetic energy per nucleon. 



Ratio 


Events 


Fit Results 




SA/A 


ACorr 


FCorr 


^He/^He 


18,035 


0.174 ±0.009 


67.7/34 


12.7% 


1.02 


0.98 


^Li/^Li 


1,046 


0.951 ±0.086 


16.1/23 


13.6% 


0.97 


0.99 


7Be/9+iOBe 


400 


1.512 ± 0.238 


19.9/23 


13.8% 


0.96 


0.99 


lOB/iiB 


1,598 


0.494 ± 0.060 


28.0/29 


13.9% 


0.96 


0.99 



tion in the Z ~ 4 sample from adjacent charges that led to 
larger systematic errors in the beryllium measurement. This 
channel was also limited by the inability to separate ^Be from 
^°Be. While for the other ratios it can be safely assumed that 
each charged species is composed of only two long-lived iso- 
topes, a few percent of ^"Be has been measured in the cosmic 
ray flux in addition to the more abundant isotopes ^Be and 
^Be (Hams et al. 2004; Webber et al. 2002). We therefore de- 
termined the ratio 'Be/(^Be-i-^'^Be) simultaneously with the 
additional parameter ^"^Be/^Be in our composition fit . The 
latter was accounted for the proper determination of the ra- 
tio ^Be/(^Be-i-^°Be) and its corresponding error. As shown in 
Fig. 8, the ratio ^°Be/^Be is poorly constrained by the data 
(between ~0 and ^^0.6 within l-tr of uncertainty). Fig. 8 
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TABLE 5 

Results for the isotopic ratios and fluxes at the 
top of instrument. 



Energy" 


J. IC/ i. J.^ IVCILILJ 


He Flux'' 


■*He Flux 


0.20 - 


0.30 


0.137 ± 0.014 


23.3 ± 2.6 


170 ± 19 


0.30 - 


0.44 


0.163 ± 0.011 


24.7 ± 2.2 


152 ± 14 


0.44 - 


0.64 


0.178 ± 0.011 


21.6 ± 1.8 


121 ± 10 


0.64 - 


0.95 


0.203 ± 0.012 


17.5 ± 1.4 


86.5 ± 6.9 


0.95 - 


1.40 


0.215 ± 0.014 


12.2 ± 1.0 


56.6 ±4.7 


Energy" 


^H/^'He Ratio 


Ratios in 0.2-1.4 GeV/n 


0.20 - 


0.30 


0.183 ±0.024 


3He/*He 


0.173 ±0.010 


0.30 - 


0.44 


0.190 ±0.020 


^LiALi 


0.912 ± 0.090 


0.44 - 


0.64 


0.188 ±0.021 


7Be/9+iOBe 


1.450 ±0.242 


0.64 - 


0.95 


0.204 ± 0.027 


iOb/Hb 


0.4695 ± 0.061 



" Kinetic energy is given units of GeV nucleon ^ . 
Fluxes are given in units of nucleon/GeV/s/m^/sr. 



also shows that the uncertainty in the ^"Be/^Be ratio has no 
dramatic consequences in the ^Be/(^Be+^°Be) ratio, given 
the weak correlation of the two parameters. The contribution 
from inelastic collisions and fragmentation was estimated as 
described in §3.2.2. Similar values as for heUum ('^ 2—3%) 
were found for both the effects. Finally, the errors from the 
TOF energy resolution and from the MC acceptance estima- 
tion were smaller than 1%. The total error assigned to the 
measurements was obtained by the sum in quadrature of all 
the noted contributions. A detailed summary is provided in 
Table 3. 

4. RESULTS AND DISCUSSION 

In the previous sections, we have described the analy- 
sis procedure adopted for the determination of the ratios 
^He/^He, ^Li/^Li, ^Bd^Be+^'^Be) and i"B/"B. The TOI 
corrections turned out to be of the same order of magni- 
tude as the estimated uncertainties, hence, the gross features 
of the measured ratios were apparent directly from the fits 
on the mass distributions. The error from the fitting pro- 
cedure was considerably larger than the other contributions. 
The most important limitations were the mass resolution (for 
He) and the limited statistics (for Li-Be-B). In particular, the 
mass resolution was limited by multiple scattering (affect- 
ing SR/R at ^0.2 GeV nucleon^^) and the TOF resolu- 
tion (affecting (5/3//3 at GeV nucleon^^). 
The results with all corrections applied are presented in Ta- 
ble 5. Results for the isotopic ratio ■^He/'*He as a func- 
tion of the kinetic energy per nucleon are shown in Fig. 9 
between 0.2 and 1.4 GeV nucleon^^ (filled circles). The 
error bars represent the total errors as discussed in §3.3. 
The figure also shows the existing data between 0.1 and 10 
GeV nucleon^ ^ measured by the balloon borne experiments 
BESS (Wangetal. 2002), IMAX (Reimeretal. 1998), the 
first flight of SMILI (Beatty et al. 1993), Hatano et al. (1995) 
and Webber & Yushak (1983). Among these, our data are the 
only data collected directly in space. Our results agree well 
with data collected by BESS in its first flight in 1993. 

Fig. 10 shows our results for Li-Be-B. The AMS-01 
data are compared with measurements made by the space 
experiments CRIS on ACE (de Nolfo et al. 2001), VOY- 
AGER (Webber et al. 2002), ULYSSES (Connell 1998), 
ISEE 3 (Wiedenbeck & Greiner 1980), and with balloon data 
from ISOMAX (Hams et al. 2004), IMP 7/8 (Garcia-Munoz 
1977) and GSFC (Hagenetal. 1977). Our results are con- 



sistent with these data within uncertainties, in particular with 
ISOMAX. 

In Fig. 11a we report the ''He and "'He differential spec- 
tra. These spectra are obtained by the combination of the 
■'He/He and He/He fractions, directly given by the •'He/''He 
ratio, with the AMS-01 helium spectrum previously pub- 
lished in Alcaraz et al. (2000). The ''He and "'He data points 
and their errors were extracted through a logarithmic inter- 
polation of helium data along our energy points. An addi- 
tional 1 % of error was added due to the interpolation proce- 
dure. Similarly, the resulting "'He spectrum has been further 
combined with the galactic deuteron spectrum published in 
Aguilar et al. (2002), to extract the ratio between deuterons, 
^H, and their main progenitors "'He. The AMS-01 data 
analysis of ^H is described in §4.6 of Aguilar et al. (2002), 
where the extraction of the deuteron signal from the vast pro- 
ton background is quantitatively discussed and the absolute 
deuteron spectrum is presented in different geomagnetic lati- 
tude ranges. The resulting ^H/^He ratio is shown in Fig. 1 lb 
together with the previous experiments BESS (Wang et al. 
2002), IMAX (de Nolfo et al. 2000), and Webber & Yushak 
(1983). While all the measurements give larger ^H/''He ratios 
than the model predictions (see below) by up to a factor 2, our 
results again show good agreement with the data from BESS. 
These data are also reported in Table 5. 

To describe our data, we show in all plots the model 
calculations of the conventional reacceleration model used 
in standard methodologies and extensively described else- 
where (Strong & Moskalenko 1998). Calculations have been 
made with the package GALPROP-v50 . 1^. 

GALPROP solves the diffusion-transport equation for a 
given source distribution and boundary conditions for the 
galactic CRs, providing steady-state solutions for the local in- 
terstellar spectra (LIS) for all the charged CRs up to Z = 28. 
The diffusion of cosmic rays through the magnetic halo is 
described by means of a rigidity-dependent diffusion coef- 
ficient D = 13 Dq {R/RqY , where Dq and Rq fix the nor- 
malization, and the spectral index S drives its rigidity depen- 
dence. The reacceleration of charged particles due to scat- 
tering on hydromagnetic waves is described as a diffusion in 
momentum space. This process is controlled by the Alfven 
speed of plasma waves moving in the interstellar medium, 
va- The code also describes energy losses due to ionization or 
Coulomb scattering, and catastrophic losses over the galactic 
disk, making use of a large compilation of cross section data 
and decay rates. To decouple all the transport equations, the 
fragmentation network starts with the heaviest nucleus and 
works downward in mass, processing primary and all sec- 
ondary nuclei produced by the cascade. This loop is repeated 
twice. 

In the parameter setting considered here, no tuning was 
done to our isotopic data. The nucleon injection spectrum 
is taken as a "broken" power law in rigidity to better match 
our total helium and proton spectra from Aguilar et al. (2002). 
Two indices I'l and 1^2 were used below and above Rb- The 
cross section database was extended using the updated cross 
section list from the version v54 (Vladimirov et al. 2011), 
which includes the production of ^H and ^He from fragmen- 
tation of heavier isotopes. The transport parameters Dq, 6 
and VA are consistent with our B/C ratio from Aguilar et al. 
(2010). In our description, we used a cylindrically sym- 

^ http://galprop.stanford.edu 
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Fig. 9. — Results for the ratio ^He/'^He between 0.2 and 1.4 GeV nucleon ^ of kinetic energy per nucleon. Other data are from BESS (Wang et al. 2002), 
IMAX (Reimer et al. 1998), Hatano et al. (1995), SMILI (Beatty et al. 1993), Webber & Yushak (1983). The dashed hne is the model calculation for the LIS 
ratio obtained with GALPROP (Strong & Moskalenko 1998). The modification of this by solar modulation for 1998 June is shown by the solid line. 
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Fig. 10. — Results for the ratios '^LiFU, ^Be/('^Be+i"Be) and ^'^B/^^B between 0.2 and 1.4 GeV nucleon"^ of kinetic energy per nucleon (•=AMS- 
01 present work, 0=CRIS/ACE (de Nolfo et al. 2001), ■=ISOMAX (Hams et al. 2004), A=VOYAGER (Webber etal. 2002), T=ULYSSES (Connell 1998), 
x=ISEE 3 (Wiedenbeck & Greiner 1980), □=IMP 7/8 (Garcia-Munoz 1977) and *=GSFC (Hagen et al. 1977)). The dashed (solid) lines are the model calcula- 
tions for the LIS (solar modulated) ratios obtained with GALPROP (Strong & Moskalenko 1998). 



metric model of the galactic halo with radius Ti=30 kpc and height Zh= 4 kpc. The relevant parameters are reported 
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Fig. 1 1 . — a) Differential spectra of '^He and ''He. Our data liave been derived from our earlier published helium spectrum combined with the "^He/^He ratio 
of this work. EiTors ai'e smaller than the size of the circles, b) The ratio ^H/^He derived from our earlier work on flux. In both figures, other data are from 
balloon borne experiments (Wang et al. 2002; Webber & Yushak 1983; Webber et al. 1991 ; Leech & O'Gallagher 1978; de Nolfo et al. 2000). The dashed (solid) 
lines are the model calculations for the LIS (solar modulated) spectra obtained with GALPROP (Strong & Moskalenko 1998). 



TABLE 6 
Propagation parameter set. 



Parameter 


Name 


Value 


Injection, break value 


Rb [GV] 


9 


Injection, index below Rg 




1.80 


Injection, index above R b 


vi 


2.35 


Diffusion, magnitude 


Do [cm^ s-i] 


5 ■ lO^s 


Diffusion, index 


(5 


0.41 


Diffusion, ref. rigidity 


flo [GV] 


4 


Reacceleration, Alfven speed 


VA [km s^^] 


32 


Galactic halo, radius 


7^ [kpc] 


30 


Galactic halo, height 


zh [kpc] 


4 


Solar modulation parameter 


4> [MV] 


450 



in Table 6; the remaining specifications are as in the file 
galdef _5 0p_5 9 9278 provided with the package. 

In the ratios of Fig. 9, 10 and 11, local interstellar (dashed 
lines) and heliospheric propagated (solid lines) calculations 
are shown. The heliospheric modulation is treated using 
the force field approximation (Gleeson & Axford 1968), with 
(/) = 450 MV as the modulation parameter to characterize the 
modulation strength for 1998 June. This is also in accordance 
with the study performed in Wiedenbeck et al. (2009) over the 
full solar cycle 23. It should be noted, however, that the force 
field approximation has no predictive power, i.e., the value 
employed for the parameter is contextual to the propagation 
framework adopted to predict the interstellar spectra of the 
CR elements. For instance, the deuteron flux in Aguilar et al. 
(2002) was described using — 650 MV after assuming a 
pure power-law energy spectrum in the LIS. 

Though large uncertainties are still present in the helio- 
spheric propagation, the general trend is that higher modu- 
lation levels correspond to lower values of the ^He/^He ra- 
tio. AMS-01 and BESS data come from periods of rela- 
tively quiet solar activity as do data from Webber & Yushak 
(0 w 400 - 650 MV). In particular, the periods of 1998 June 
(AMS-01 flight) and 1993 July (BESS flight) were character- 
ized by very similar solar conditions according to the sunspot 
data (Temmer et al. 2002) and can be directly compared. 
Stronger modulations were present when IMAX (0 « 700- 



- 850 MV) and SMILI {(i> w 1200 - 1300 MV) were active. 

In summary, the secondary to primary ratio '^He/^'He, which 
is much more sensitive to the propagation parameters, seems 
to be well described by the model under these astrophysical 
assumptions. The model is also consistent with the Z > 2 
ratios of Fig. 10, though large errors are present in the current 
data. As the Li-Be-B elements are of secondary origin, these 
ratios are less sensitive to the galactic transport, and may be 
useful to investigate the nuclear aspects of the CR propagation 
(fragmentation, decay, and breakup). On the contrary, our 
^H/ He data of Fig. lib give a larger ratio than the model 
predictions, and this tendency is also apparent from the other 
experiments. Understanding this possible discrepancy may 
require a thorough investigation of the production cross 
sections, in particular for the reactions induced by cosmic ray 
protons and helium nuclei. 

5. CONCLUSIONS 

The AMS-01 detector measured charged cosmic rays dur- 
ing 10 days aboard the Space Shuttle Discovery in 1998 June. 
Owing to the large number of helium events collected and the 
absence of the atmospheric effects, we have precisely deter- 
mined the ratio ''He/''He in the kinetic energy range from 0.2 
to 1.4 GeV nucleon^^. The average isotopic ratios ^Li/^Li, 
^Be/(^Be-i-^''Be) and ^"^B/^^B have been measured in the same 
energy range. The ratio -^H/^He and the spectra of ■^He and 
^He are also reported. Our results agree well with the previous 
data from BESS and ISOMAX and can provide further con- 
straints to the astrophysical parameters of cosmic ray propa- 
gation. In the analysis procedure adopted in this work, Monte 
Carlo simulations were essential for understanding the instru- 
ment performance, its acceptance, the role of interactions and 
for modeling the mass distributions. We expect, with AMS- 
02, to achieve much more precise results over wider energy 
ranges. 
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